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O- ' ABSTRACT 

C") \ We present B, R, and /-band selected galaxy catalogues based on the Munich Near- 

Infrared Cluster Survey (MUNICS) which, together with the previously used K- 
selected sample, serve as an important probe of galaxy evolution in the redshift range 
< z < 2. Furthermore, used in comparison they are ideally suited to study selection 
■ effects in extragalactic astronomy. The construction of the B, R, and /-selected pho- 

Q^' tometric catalogues, containing ~ 9000, ~ 9000, and ~ 6000 galaxies, respectively, is 

described in detail. The catalogues reach 50% completeness limits for point sources 
£5 ' of B ~ 24.5 mag, R ~ 23.5 mag, and / ~ 22.5 mag and cover an area of about 0.3 

^ \ square degrees. Photometric redshifts are derived for all galaxies with an accuracy of 

C3 • + z ) — 0.057, very similar to the Zsf-selected sample. Galaxy number counts in 

the B, V, R, I, J, and K bands demonstrate the quality of the dataset. 

The rest-frame colour distributions of galaxies at different selection bands and 
. redshifts suggest that the most massive galaxies have formed the bulk of their stellar 

$— i ' population at earlier times and are essentially in place at redshift unity. We investigate 

the influence of selection band and environment on the specific star formation rate 
(SSFR). We find that ZT-band selection indeed comes close to selection in stellar 
mass, while Z3-band selection purely selects galaxies in star formation rate. We use a 
galaxy group catalogue constructed on the ZT-band selected MUNICS sample to study 
possible differences of the SSFR between the field and the group environment, finding 
a marginally lower average SSFR in groups as compared to the field, especially at 
lower redshifts. 

The field-galaxy luminosity function in the B and R band as derived from the 
/?-selected sample evolves out to z ~ 2 in the sense that the characteristic luminosity 
increases but the number density decreases. This effect is smaller at longer rest-frame 
wavelengths and gets more pronounced at shorter wavelengths. Parametrising the 
redshift evolution of the Schechter parameters as M*(z) = M*(0) + a ln(l + z) and 
$*(z) = <1>*(0)(1 + z) b we find evolutionary parameters a ~ —2.1 and b ~ —2.5 for 
the B band, and a ~ — 1.4 and b a — 1.8 for the R band. 

Key words: surveys - galaxies: evolution - galaxies: fundamental parameters - 
galaxies: luminosity function - galaxies: photometry - galaxies: stellar content 
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1 INTRODUCTION 

Investigating the evolution of galaxies with time requires the 
analysis of statistical properties of galaxies at various cos- 
mic epochs. Studies of the luminosity function of galaxies at 
different wavelengths, of the galaxies' mass function or star 
formation rates at different redshifts are important meth- 
ods to address the problem of the formation and evolution 
of galaxies within observational cosmology. These investiga- 
tions rely on galaxy surveys large enough to yield sufficient 
statistical accuracy and deep enough to trace the properties 
of the galaxy population to high redshifts and faint intrinsic 
luminosities. 

Since different wavelengths sample light from different 
stellar populations within the galaxies, the selection band 
of a galaxy survey is of fundamental importance. At the 
reddest wavelengths, a galaxy's light is dominated by old 
stellar populations and hence reflects its total stellar mass 
l|Rix fc Riekell993r ). Moving to bluer wavelengths, light from 
younger stellar populations becomes important, and samples 
selected in bluer filters become biased toward on-going star 
formation activity. 

Historically, by the 1980s the increased number of 4-m 
class telescopes and advances in electronic detector technol- 
ogy enabled deep surveys of the universe out to redshifts 
z ~ 1, when the universe had on ly about half of its present 
age (see IColless|[l997l ; lEllisll 19971 for reviews). 

Traditionally, many surveys were conducted in the B 
band using photographic plates at Schmidt telescopes. When 
multi-object spectrographs became available at 4-m class 
telescopes, many groups started follow-up observations of 
these surveys. One noteworthy redshift survey of that kind 
is the Autofib Redshift Survey probing the evolution of the 
luminosity function o ut to redshifts z ~ 0.8 with roug hly 
1700 galaxy redshifts (|Ellis et al.lll996l ; iHevl et al.lll997T l. 

In the middle of the 1990s, the I-band selected Canada- 
France Redshift Survey (CFRS) marked an important point 
in the study of galaxy evolution with redshift surveys. The 
selection in the I band made it possible to study the evolu- 
tion of the evolved, massive galaxies rather than star forming 
galaxies picked up in the more traditional B-band selected 
samples. The CFRS comprises spectra of more than 1000 ob- 
jects with 17.5 ^ Iab ^ 22.5, 591 of which are galaxies with 
secure redshifts in the range ^ z ^ 1.3 and en abled stud- 
ies o f the evolution of the luminosity function l|Lillv et al.l 
1 19951 ) and of the lumino sity density and s tar formation rate 
density of the universe (|Lillv et al.1 1 1996ft . finding very lit- 
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Table 1. The ten MUNICS fields. The table gives the field name, 
the field coordinates (right ascension a and declination 5 as well 
as galactic coordinates I and b) for the equinox 2000.0 and the 
effective area of the field in square arc-minutes. The effective area 
is the area of the sky covered by observations in all six filters. 



Field 


a 


6 


I 


b 


Area 




(2000.0) 


(2000.0) 


[arcmin 2 ] 


S2F1 


03:06:41 


+00:01:12 


178.66 


-47.67 


119.1 


S2F5 


03:06:41 


-00:13:30 


178.93 


-47.84 


123.7 


S3F1 


09:04:38 


+30:02:56 


195.37 


+40.64 


116.6 


S3F5 


09:03:44 


+30:02:56 


195.32 


+40.45 


107.6 


S4F1 


03:15:00 


+00:07:41 


180.58 


-46.05 


105.6 


S5F1 


10:24:01 


+39:46:37 


180.85 


+57.04 


115.7 


S5F5 


10:25:14 


+39:46:37 


180.76 


+57.27 


105.3 


S6F1 


11:55:58 


+65:35:55 


131.90 


+50.55 


117.1 


S6F5 


11:57:56 


+65:35:55 


131.60 


+50.62 


116.5 


S7F5 


13:34:44 


+16:51:44 


349.45 


+75.66 


119.1 


All 










1146.3 



tie evolution in the luminosity and number density of red 
galaxies in the redshift range < z < 1. 

Selection in near-infrared filters like the _K"-band comes 
even closer to a selection in stellar mass. When suitable near- 
infrared detectors became available at large telescopes, a 
number of K-h&nd selected surveys were undertaken which 
can probe the evolution of massive galaxies out to redshift 
z ~ 1. Examples for s uch surveys are the Hawaii Deep 
Field s (|Cowie et al.lll994l ). the K20 Survey l|Cimatti et all 
12002ft and the Mun ich Near-Infrared Cluster Survey (MU- 
NICS; [DroaelaD 2001b). 

In this paper we present B, R, and I band selected 
MUNICS catalogues which can be used for comparison with 
previous work and, together with the K-band selected cat- 
alogue, to study selection effects in extra-galactic surveys. 
They are an important tool to trace different parts of the 
field galaxy population, and to discriminate between selec- 
tion biases and real evolutionary effec ts. Furtherm ore, this 
paper also serves as a general update to lDrorv et al] (2001b) . 
Since the publication of that paper, the survey's high-quality 
multi-colour imaging part has not only grown in area, but 
B-band imaging has become available, and the catalogue 
has been improved in a number of ways. 

This paper is organised as follows. Section [2] gives a 
brief overview of the Munich Near-Infrared Cluster Survey 
(MUNICS), while Section [3] describes the construction and 
content of the optically selected catalogues presented in this 
paper. In Section [4] we present galaxy number counts in all 
six filters and compare them to previous studies. Section [5] 
discusses results on the rest-frame colour distributions of 
galaxies in the different catalogues, Sections [6] and [7] elabo- 
rate on the influence of selection effects and environment on 
the SSFR, respectively, while Section [8] presents results on 
the evolution of luminosity functions, before we summarise 
our findings in Section [9] Throughout this work we assume 
Q M = 0.3, Qa = 0.7 and H = 70kms~ 1 Mpc -1 . All mag- 
nitudes are given in the Vega system. 
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Figure 1. Result of Monte-Carlo completeness simulations with artificial point sources placed in the B (upper left-hand panel), R (upper 
right-hand panel), I (lower left-hand panel), and i^-band images (lower right-hand panel) of MUNICS. 



2 THE MUNICH NEAR-INFRARED CLUSTER 
SURVEY (MUNICS) 

The Munich Near-Infrared Cluster Survey (or MUNICS for 
short) is a wide-field medium-deep imaging s urvey in the 
near-infrared and optical initially described in iDrorv et al.l 
(2001b). Dedica ted follow-up spect roscopy is available for 
- 600 galaxies l|Feulner et al.ll2003h . The main part of the 
survey consists of 10 fields the details of which are sum- 
marised in Table [1] For all these fields photometry in K' , 
J, I, R, V, and B is available, with limiting magnitudes 
ranging from K' ~ 19.5 to B ~ 2 4.5 (50% completeness 
for point sources; ISnigula et aLl l2002). The total area of this 
part of the survey is 1146.3 square arcmin (or 0.32 square 
degrees). However, the fields S3F1 and S4F1 have less than 
satisfactory quality, especially concerning their photomet- 
ric calibration. Since an accurate calibration of photometric 
colours is essential for high-quality photometric redshifts, 
we exclude these two fields from any further analysis (as 



has been done with the if-selected sample). The total area 
of the eight high-quality fields is then 924.1 square arcmin 
(or 0.26 square degrees). 



Most of the research on field galaxy evolution and 
galaxy clusters in the MUNICS project has bee n carried out 
with the Jf-band selected catalogue presented in lDrorv et al.l 
(2001b) . This includes studies of luminosity f unction evo- 
lution (|Feulner et all 120031; [Drorv et all 120031), th e stellar 
mass function of galaxies ( Drorv et al.l 2001a, |2004f). the in 



tegra ted specific star formation rate of groups (Fculncr et al 



2006), and the cluste r catalogue described in lBotzler et al 



l|2007ft . lFeumer et all (2005b), however, a study of the evolu- 
tion of the specific star- formation rate with redshift, is based 
on the /-band selected catalogue. In this paper, the R-, 
and _B-selected MUNICS samples will be described for the 
first time in detail. 
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Early types Late types 




^ [H™] SED model number 

Figure 2. Left-hand panel: Template spectral energy distributions (SEDs) for galaxies used in the computation of photometric redshifts 
for the /-, R-, and B-selected samples. Right-hand panel: The distribution of SED types in the four MUNICS catalogues. The 
SED number is an internal SED identification: Small numbers refer to redder (early-type) galaxies, higher numbers to bluer (late-type) 
galaxies. 



3 THE CONSTRUCTION OF OPTICALLY 
SELECTED CATALOGUES 

In this section we want to discuss the construction and the 
properties of the B, R, and /-band selected photometric cat- 
alogues as well as the measurement and reliability of pho- 
tometric redshifts. We will often refer to the different cat- 
alogues as MUNICS_B, MUNICS_R, MUNICSJ, and MU- 
NICS_K for short. 



3.1 Object Detection 

As for the ^-selected catalogue, object detection on the B, 
R, and /-band image s was perform ed using the YODA source 
extraction software l|Drorv| feoo&t ) . Sources are detected by 
requiring a minimum number N pix of consecutive pixels to 
lie above a certain threshold t expressed in units of the local 
RMS a of the background noise. To ensure secure detection 
of faint sources, the images are convolved with a Gaussian of 
full width at half maximum (FWHM) similar to the seeing 
in the image. The choice of the number of consecutive pixels 
Npi x and the threshold t is a compromise between limiting 
magnitude at some completeness fraction, say 50 per cent, 
and the nu mber of to lerable spurious detecti ons per unit 
image area (|Sahalll995l . see also the discussion in lDrorv et al.l 
2001b). 

To find reasonable values for N P i x and t we performed 
simulations with varying minimum number of consecutive 
pixels on the optical images of one of our mosaic fields 
(S6F5). Instead of fixing two parameters (the threshold and 
the convolution kernel) and keeping only one free parame- 
ter (as with the Zf-band detection), we chose only one fixed 
parameter: The size of the convolution kernel was set to be 
equal to the (Gaussian) point-spread function (PSF) of the 
images. Thus we kept both the threshold and the number 
of consecutive pixels as free parameters in our simulations. 



Choosing a PSF-like Gaussian kernel for the image convolu- 
tion is necessary to ensure secure detection of faint, compact 
objects in the images. The number of false detections was 
measured by running the detection algorithm on the same 
images multiplied by —1. 

From these simulations, the best choice of parame- 
ters is t — 4a for the detection threshold, and N P i X = 
-k(FWHM/2) 2 for the minimum number of consecutive pix- 
els required for an object. With these settings, the contami- 
nation rate (the rate of false detections) is actually below 1 
per cent. 

Since one of the ideas of this paper is a comparison of 
properties of the /('-band and the optically selected galaxies, 
we have convinced ourselves that this change of detection pa- 
rameters does not influence the results presented here. The 
reason is that different detection parameter choices mostly 
affect the object statistics at the very faintest apparent mag- 
nitudes. Since these objects also have large photometric er- 
rors, they are usually excluded from any further analysis and 
thus cannot influence the results of our work. 

In Figure Q] we show the completeness functions for the 
B, R, and /-band images of the eight high-quality MUNICS 
mosaic fields, based on Monte-Carlo simulations of artificial 
point-sources placed at random in the images. In these simu- 
lations, 200 artificial objects with a Moffat-type PSF having 
the same FWHM as stars in the frames and a constant mag- 
nitude distribution were added to the images using the iraf 
ARTDATA package. After running the object detection algo- 
rithm on the images applying the same detection parameters 
as for the main catalogue, the fraction of re-detected arti- 
ficial objects as a function of magnitude is computed. This 
procedure was repeated 500 times in order to decrease sta- 
tistical errors. An extensive discussion of completeness sim- 
ulations for ext ended objects in MU NICS A'-band images 
can be found in lSnigula etaH (I2002T I. 

One can see from the diagram that the 50% complete- 
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Figure 3. Upper panels: A comparison of photometric and spectroscopic redshifts for MUNICSJ3 (left), MUNICS_R (middle), and 
MUNICSJ (right). Filled symbols are galaxies, open symbols quasars. Lower panels: The corresponding histograms of the rcdshift 
differences toget her with a Gaussia n fit showing the quality of the photometric redshifts. The corresponding diagrams for MUNICS_K 
are published in lDrorv et al. (2003). 



ness limits for point sources as derived from these simula- 
tions are B ~ 24.5 mag, R ~ 23.5 mag, I ~ 22.5 mag, and 
K ~ 19.5 mag. 



3.2 Photometry 

Photometry was done in elliptical apertures the shape of 
which was determined from the first and second moments 
of the light dist ribution in the detection image, as described 
in iDrorv (|2003l ). and additionally in fixed size circular aper- 
tures of 5 and 7 arc seconds diameter. To ensure measure- 
ment at equal physical scales in every pass-band, the in- 
dividual frames were convolved to the same seeing FWHM, 
namely that of the image with the worst seeing in each field. 

Aperture fluxes and magnitudes were computed for each 
object present in the B, R, and /-band catalogues irrespec- 
tive of a detection in any other band. For this purpose the 
centroid coordinates of the sources found in the detection 
images were transformed to the other frames using geomet- 
ric transformations between the image coordinate systems. 
For this purpose, the images in the other five filters of each 
field were registered against the image in the detection band 
by matching the positions of ~ 200 bright homogeneously 
distributed objects in the frames and determining the coor- 
dinate transform from the detection system to each image 
in the other five pass-bands using the tasks xyxymatch 
and GEOMAP within iraf. The scatter in the determined so- 



lutions is less than 0.1 pixels RMS in the transformation 
from / or B to the other optical bands, and less than 0.2 
pixels RMS from I or B to the near-infrared frames. Note 
that the frames themselves are not transformed to avoid 
artefacts from the re-sampling. We only determine accurate 
transformations and apply these later to the apertures in 
the photometry process, where the shapes of the apertures 
were transformed using only the linear terms of the trans- 
formation. 

As a result of this method, small differences in object 
centring and differences in the object's shape in the differ- 
ent detection images lead to slightly different photometry in 
the various catalogues. Nevertheless, a direct comparison of 
magnitudes in the same filter, but derived in the different 
selection catalogues is an important consistency check. The 
comparison shows excellent agreement between MUNICS_B, 
MUNICS _R, M UNICSJ, and M UNICS JK (comparison plots 
are published in iFeulneijliool) . 

3.3 Photometric Redshifts 

Photomet ric redshifts are der ived using the method pre- 
sented in iBender et alj (2001), a template matching algo- 
rithm rooted in Bay esian st a tistics closely resembling the 
method presented bv lBem'ted (|200(t ). The application of this 
method to the -selected MUNI CS sample is described in 
great detail in lDrorv et~aT] (|2003f ). 
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Z phot Z phot 

Figure 4. Plots of the absolute magnitudes Ms, Mr, Mj, and Mk versus the photometric redshift z p hot for MUNICS-B [upper left-hand 
panel), MUNICS-R {upper right-hand panel), MUNICSJ {lower left-hand panel) and MUNICS-K {lower right-hand panel), respectively. 
The different colours denote the different model SEDs, where the colour corresponds roughly to the model colour, i.e. early-type galaxies 
are shown in red, late-type galaxies in blue. Objects spectroscopically classified as active galactic nuclei (AGN) are additionally marked 
as open squares. The dashed lines represent the expected absolute magnitudes as a function of redshift for an early-type SED {red dashed 
line) and a very late type SED {blue dashed line) at the limiting magnitudes of the samples, i.e. B = 24.5, R = 23.5, / = 22.5, and 
K' = 19.5. 



The left-hand panel of Figure [2] shows the final tem- 
plate spectral energy distribution (SED) library used to de- 
rive photometric redshifts in what follows. It is the same 
set of SEDs also used for the K-band selected catalogue, 
since these SEDs already cover a vast range from very red 
old models to young star-bursting models. As we show in 
the right-hand panel of Figure [2] the difference between the 
differently selected catalogue lies then in the distribution 
of selected SEDs: While in the if-selected sample the algo- 
rithm picks preferentially red galaxy types (and only very 
few heavily star- forming objects), the distribution for the 
/- and /^-selected samples is more balanced, and indeed re- 
versed for the _B-selected sample. 

Figure [3] compares photometric and spectroscopic red- 
shifts for all ~ 600 objects with spectroscopic redshifts 
within the MUNICS fields and shows the distribution of red- 
shift errors. The typical scatter in the relative redshift error 



Az/{1 + z) is 0.057. This is very similar to the results for th e 
K-baiui selected catalogue described in iDrorv et all (|2003h . 
The mean redshift bias is negligible. The distribution of the 
errors is roughly Gaussian. There is no visible difference be- 
tween the distributions among the survey fields. Although 
this performance is encouraging, it is important to say that 
the spectroscopic data become sparse at z > 0.6 and there 
are only very few spectroscopic redshifts at z > 1. Note that 
many quasars are among the few dramatic outliers which 
is to be expected from the power-law like SEDs of active 
galactic nuclei. 

Many objects in the MUNICS fields are, of course, de- 
tected in more than one band. As a consistency check, we 
present comparisons of the photometric redshifts of these ob- 
jects derived in different catalogues in Appendix[X] showing 
excellent agreement. 

In Figure [4] we show the distribution of absolute B, 
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Figure 5. Redshiftcd spectral energy distributions (arbitrarily normalised) of an early-type galaxy (thick line) and a star-burst galaxy 
(thin line) and the selec tion filters B, R, I, and K (shaded areas) for various redshifts. The early-type SED is constructed from 
iBruzual &i Charlo t (2003) models and has an age comparable to the age of the universe at each redshift. 
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Figure 6. Histograms of photometric redshifts for MUNICSJ3 (dotted line), MUNICS.R (dashed line), MUNICSJ (dot-dashed line), 
and MUNICS JC (solid line). 



R, I, and K magnitude versus photometric redshift for the 
four MUNICS samples and for different model SEDs, rang- 
ing from early types (redder colours) to late types (bluer 
colours). We also show the expected curve for a early-type 
model SED (red dashed line) and a very late type model 
SED (blue dashed line) at the limiting magnitudes of the 
samples, i.e. B = 24.5, R = 23.5, / = 22.5, and K' = 19.5. 

The object distribution is generally in very good agree- 
ment with the expectations. Outlying objects with suspi- 
ciously high luminosities are likely to be partly photometric- 
redshift outliers and active galactic nuclei (AGN). Indeed, 
spectroscopically identified AGN are marked in the Figure, 
and they are found in this region of the diagram. Note that 
the spectroscopic follow- up of MUNICS wa s pre-selected 
against point sources (see lFeulner et al.1l2003l ). and that the 
existing AGN spectra are mostly drawn from a sample of 



AGN selected for their red colour (Wisotzki et al., in prepa- 
ration). The lack of objects with z ^ 1.5 in MUNICSJ seems 
puzzling at first. From their apparent magnitudes we should 
be able to detect blue galaxies at these redshifts. To inves- 
tigate this apparent failure, we have selected objects in this 
redshift range from MUNICS Jl and plotted their position in 
the /-band images. It turns out that these objects are simply 
not detected due to the higher noise level in the /-band im- 
ages caused by fringing. The higher noise level reduces the 
overall depth of the sample and affects galaxies at all red- 
shift. Looking at this effect in terms of the luminosity func- 
tion (LF) at different redshifts, the incompleteness cuts away 
the faint end of the LF in each redshift bin, with the limit 
moving to higher luminosities with increasing redshift. At 
redshift z ~ 1.5 the /-band magnitude limit of MUNICSJ 
reaches the bright end of the LF, i.e. only very few galax- 
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ies in this redshift range are detected. We have verified this 
by taking the /^-selected catalogue and plotting that /-band 
magnitude histogram of all galaxies beyond z = 1.5 which 
indeed falls below the /-band detec tion limit. Also, th e much 
deeper /-selected FDF catalogue jHeidt et alj|2003h cut at 
the MUNICS J limit shows a redshift distribution similar to 
the one for MUNICSJ. Note that this lack of galaxies be- 
yond z ~ 1.5 in the /-selected sample does not cause any 
problems as long as one excludes these higher redshifts from 
any analysis. 

The distribution of SED types in Figure [4] clearly il- 
lustrates the importance of selection effects at higher red- 
shift: While the /^-selected sample traces the population of 
early-type galaxies out to redshifts beyond z ~ 1.5 (at the 
limiting magnitudes of MUNICS; indeed, that is how the 
survey was designed), it lacks late-type galaxies at all red- 
shifts. The /^-selected sample on the other hand does not 
contain early-type galaxies beyond z ~ 0.8 (again at the 
limiting magnitudes of MUNICS), but is much better suited 
for tracing late-type galaxies at low and high redshift. 

These selection effects can be easily understood by look- 
ing at the position of redshifted SEDs with respect to the 
four detection filters as shown in Figure [5] where we plot 
redshifted SEDs for an early-type and a star-burst galaxy. 
Early-type galaxies disappear from a /3-selected catalogue 
at z ~ 1 and from R as well as /-selected catalogues at 
2 ~ 1.5, while they remain visible to much higher redshifts 
in MUNICSJK. Similarly, the detection of very blue galax- 
ies at redshifts z ~ 3 in MUNICS_B and MUNICS_R can be 
explained by the strong rest-frame ultraviolet emission at 
wavelengths longer than that of the Lyman break entering 
the detection bands. 

The photometric redshift histograms for MUNICS_B, 
MUNICS_R, MUNICSJ, and MUNICSJK can be found in 
Figure [6] Note the high redshift tails in the distributions 
for MUNICSJ3, MUNICS_R, and MUNICSJK, caused by 
luminous red galaxies at 1 < z < 2 and the ultravio- 
let emission of blue galaxies shifted into the optical bands, 
respectively. 

3.4 Star— Galaxy Separation 

For computing statistical properties of the field-galaxy pop- 
ulation like its luminosity function, it is necessary to remove 
all stars from the analysis. In a photometric catalogue, this 
can be done in two ways. Either by looking at an object's 
morphology in the image, i.e. whether it appears to be sim- 
ilar to the point-spread function (PSF) of the image or ex- 
tended, or by looking at the spectral energy distribution 
(SED), in this case as traced by the six- filter photometry, 
and comparing it to template spectra of stars and galaxies. 
The disadvantage of the morphological approach is clearly 
its tendency to fail at fainter magnitudes (because of the 
lower signal-to-noise ratio and because distant galaxies look 
more and more compact), whi le it works re asonably well at 
brighter magnitudes (see, e.g.. iDrorv et alj 2001b). 

For this investigation, we decided to take the sec- 
ond route. The photometric redshift technique as described 
above compares the photometry in the six MUNICS filters 
to template SEDs of stars and galaxies. Each fit is assigned 
a % 2 value, and we can simply discriminate between stars 
and galaxies by comparing the \ 2 values of the best-fitting 




2 

X galaxy 

Figure 7. This diagram shows the loci of spectroscopically clas- 
sified galaxies (open circles) and stars (crosses) in a plane defined 
by the x 2 values of the best-fitting galactic and stellar spectral- 
energy distributions (SEDs) for MUNICS.B. The line used to 
discriminate the two classes is also indicated. The corresponding 
diagrams for MUNICS Jl, MUNICSJ, and MUNICS _K look very 
similar. 

stellar and the best-fitting galactic SED. More specifically, 
we classify all objects as stars for which Xstar < Xgaiaxy 

This procedure can be teste d by checking again st the 
spectroscopic sample described in lFeulner et al.l (j 2003). This 
is illustrated in Figure [7] where we show the \ 2 values of ob- 
jects spectroscopically classified as stars or galaxies, respec- 
tively, for the /3-selected catalogue. The diagrams for the 
other samples are very similar. Clearly, the SED classifier 
does work very well. Note that in all comparisons between 
MUNICSJB, MUNICSJ*,, MUNICSJ, and MUNICSJK pre- 
sented in this paper we have used the same x 2 -based classi- 
fication method. 



4 GALAXY NUMBER COUNTS 

Galaxy number counts, although not as widely accepted as 
important tool for the study of galaxy evolution and cos- 
mology as in the past, still serve as an interesting probe 
of a survey's data quality. Since the publicati on of galaxy 
number counts from MUNICS in IDrorv et~ai1 (2001b), the 
survey has considerably increased in size and quality. More- 
over, B-band number counts have never been published for 
the MUNICS project, so we show galaxy number counts in 
all six filters B, V, R, I, J and K in Figure [S] The counts 
were derived from detection catalogues constructed on all 
filter images. In contrast to the analysis presented later in 
this paper, star-galaxy separation is based on the morpho- 
logical classifier described in detail in IDrorv et al.l (2001b). 
The reason for this is that no photometric redshift data are 
available for the V, R and J-band detections, so we cannot 
rely on the x 2 -based classification. However, the two clas- 
sification methods agree very well. Completeness-corrected 
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galaxy number counts from MUNICS are presented here for 
the first time. The completeness correction is based on the 
simulations described in Section 13.11 All errors given are 
Poissonian errors. 

The galaxy number counts presented in Figure [8] for all 
six MUNICS filters show excellent agreement with existing 
data and clearly demonstrate the quality of our dataset. The 
completeness-corrected values for the number counts are also 
summarised in Table [2] 



5 DISTRIBUTION OF REST-FRAME 
COLOURS 

In this paper, we want to investigate the differences of the 
galaxy populations in K, I, R, and B-selected samples and 
follow their evolution with redshift. One simple, but reward- 
ing way of doing this is to study rest-frame colo ur distribu- 
tions as a function of luminosity and redshift. I Cole et all 
l|200ll ) divided their A"-selected sample of local galaxies into 
three luminosity classes and looked at the distribution of 
rest-frame bj — K s and J — K s colours. Since the A"-band 
luminosity is a measure of a galaxy's stellar mass, the lumi- 
nosity classes correspond to a binning in stellar mass. Their 
main finding is that at smaller K-luminosities (i.e. smaller 
masses) there is a more and more prominent population of 
blue, star forming galaxies. It is interesting to test this also 
for higher redshifts and different selection bands. 

The solid line in the upper left-hand panel of Fig- 
ure [5] shows the B — K rest-frame colour distributions of 
_K"-selected M UNICS gal a xies at ^ z ^ 1, confirming 
the result of ICole et all (|200ll ). Moreover, we present in 
this Figure the same distributions for /-band (dot-dashed 
line) as well as ii-band selected MUNICS galaxies (dashed 
line), showing clearly that the main difference between near- 
infrared and optically selected samples are the bluer colours 
of low-luminosity (low-mass) galaxies. This is even more pro- 
nounced for the _B-selected sample shown as a dotted line. 

Four conclusions can be immediatel y drawn from thi s 
distributions. Firstly, as already noted by I Cole et alj (|200lf ). 
a population of blue, star forming galaxies contributes more 
strongly at fainter luminosities. Secondly, this populations 
becomes more numerous at bluer selection wavelengths. 
Thirdly, the colour distribution seems to be wider at fainter 
absolute A'-band magnitudes (although part of this effect 
may be attributed to the larger errors in the rest-frame 
B — K colour for these objects). Fourthly, and maybe most 
importantly, the colour distributions of high-mass galaxies 
change very little as a function of redshift. Hence one would 
not expect much variation of the bright end of the galaxy 
luminosity function in different selection bands like K and 
7. The faint end, however, will be certainly affected in the 
sense that in bluer selection bands the faint-end slope will 
likely be steeper. 

In the other three panels of Figure[9]we present the rest- 
frame B — K colour distributions in the two redshift inter- 
vals [0.0,0.5] and [0.5, 1.0] for MUNICSJ3 (upper right-hand 
panel), MUNICS_R (lower left-hand panel) and MUNICS_K 
(lower right-hand panel). First, it is very interesting that 
the colour distributions agree so well for the highest luminos- 
ity objects (irrespective of the selection filter). Furthermore, 
there is a clear trend with redshift in the sense that lower- 



luminosity objects get bluer. That means that even at red- 
shift z ~ 1 the increased star-formation rate is dominated by 
low-luminosity (low-mass) systems. Finally, it is worth not- 
ing that this evolutionary trend with redshift is barely visible 
in the Tf-selected sample, clearly evident in the i?-selected 
sample, but gets very large in the B-selected sample: Indeed, 
for the low-mass objects with —21.75 > Mk > —23.00, 
the shift between redshift z ~ 0.25 and z ~ 0.75 is 
A(B -K)~ 0.3 for MUNICSJK, A(B - K) ~ 0.7 for MU- 
NICS_R, but as large as A(B - K) ~ 1.0 for MUNICS_B. 

How can this be interpreted? The rest-frame colour of 
a galaxy is governed by the age of its stellar population, 
its star-formation activity, the metalicity and the dust con- 
tent. Let us neglect the influence of metalicity and dust for 
the moment. The lack of redshift evolution of the rest-frame 
colours of high J\-band luminosity galaxies over < z < 1 
means that these objects must have built up their stellar 
mass at earlier times. The fact that also at lower K lumi- 
nosities galaxies of similar rest-frame colour can be found 
even in the higher redshift bin indicates that also part of 
the lower mass objects formed early, although the major- 
ity of them had (or has) major star-formation activity at 
2 < 1. It is also clear from these colour distributions that 
the rise of the star-formation rate to redshift unity is mostly 
driven by lower-mass galaxies. Note that 'lower mass' in this 
context includes intermediate-mass spirals for which there 
is evidence for in creased star-formation activity in the past 
ijBell et al.ll2005h . 

Overall, this picture of different paths of evolution for 
the most massive and the lower-mass galaxies is in good 
agreement with recent findings on the evolution of the spe- 
cific star-formation rate (SSF R), i.e. the star forma tion rate 
per unit stellar mass (see, e.g.. iFeulner et al1l2005l a and ref- 
erences therein). We will discuss the SSFR in more detail in 
Sections [6] and [7] of this paper. 



6 THE SPECIFIC STAR FORMATION RATE 
AND SELECTION EFFECTS 

In this Section we will discuss the influence of the selection 
band on the specific star formation rate (SSFR), i.e. the star 
formation rate per unit s tellar mass. 

As in IFeulner et all (2005b), we estimate the star for- 
mation rates (SFRs) of our galaxies from the SEDs by de- 
riving the luminosity at A — 2800 ± lOOA and converting 
it to an SFR as described in lMadau et al.l 1 19981) assu ming 
a Salpeter initial mass function (IMF; ISalpeterlll955T ). We 
have convinced ourselves that these photometrically derived 
SFRs are in reasonable agreement with spectroscopic indica- 
tors for objects with available spectroscopy. Note that since 
our bluest band is B, this is an extrapolation for z < 0.4. 
Hence we restrict any further analysis to redshifts z > 0.4, 
where the ultraviolet continuum at A ~ 2800A is shifted into 
or beyond the B band. 

Stellar masses are computed from the multi-colour 
photometry using a method similar to the one used in 
iDrorv et al.l ([2004). It is described in detail and tested 
against spectroscopic and dynamical mass estimates in 
Drory, Bender & Hopp l|2004h . In brief, we derive stellar 
masse s by fitting a grid of stella r population synthesis mod- 
els bv lBruzual fc Charlotl (2003) with a range of star forma- 
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Figure 8. Galaxy number counts from MUNICS (large squares) in B (upper left-hand panel), V (upper right-hand panel), R (middle 
left-hand panel), I (middle right-hand panel), J (lower left-hand panel) and K (lower right-hand panel). The filled squares indicated 
completeness-corrected counts, whereas the open squares represent the uncorrected values. Literature counts are shown for comparison. 
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Table 2. Completeness corrected galaxy number counts from MUNICS in the B, V, R, I, J and K bands, log JV and (?i os n are given, 
where N is in units of mag - Meg - 2 . 
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tion histories (SFHs), ages, metalicities and dust attenua- 
tions to the broad-band photometry. We describe star forma- 
tion histories (SFHs) by a two-component model consisting 
of a main component with a smooth SFH oc exp(— t/r) and 
a burst. We allow SFH time-scales r £ [0.1, oo] Gyr, met- 
alicities [Fe/H] 6 [—0.6,0.3], ages between 0.5 Gyr and the 
age of the universe at the objects redshift, and extinctions 
Ay 6 [0, 1.5]. The SFRs derived from this model fitting is in 
good agreement with the ones from the ultraviolet contin- 
uum. Note that we apply the extinction correction derived 
from this fitting also to the SFRs. 

Figure [TO] presents the SSFR versus stellar mass di- 
agram for four redshift bins as derived from MUNICSJK. 
The influence of the limiting apparent magnitude of the K- 
selected catalogue is clearly visible in the lower left-hand 
part of the panels. The boundary for MUNICSJK is marked 
by the red solid line, while the orange and blue lines give 
the limits of the point distributions for MUNICS J and MU- 
NICSJ3, respectively. 

The selection band clearly affects the slope of this 
boundary. While for the if-selected catalogue the line is 
steeper (i.e. closer to a selection in stellar mass), it gets grad- 
ually shallower for selection bands at shorter wavelengths. 
For MUNICSJB, the ridge runs essentially parallel to lines of 
constant SFR, demonstrating that B-band selection indeed 
selects galaxies according to their star-formation activity. 

For massive galaxies (log M» > 9.5) with non-negligible 
star- formation activity (SFR > 1 Mq yr -1 ) at redshifts 



z < 1 there is little difference between the three selection 
bands. This can also be seen in Figure QTJ where we compare 
the average SSFR as a function of redshift for galaxies in 
three different mass bins as derived from MUNICSJB and 
a combined sample from th e J-selected FDF, a nd the K- 
selected GOODS catalogue ()Feulner et all 120051 '). We chose 
MUNICSJB for this exercise because it is the only sample 
not affected by incompleteness for the lowest mass bin. Val- 
ues from MUNICS_B and FDF/GOODS-S are in remark- 
able agreement. The good agreement between the /-selected 
FDF and the if-selec ted GOODS-S catalogue has already 
been demonstrated in lFeulner et alJ (2005a). We a lso show 
local values as derived by iBrinchmann et al.l (|2004T ) for the 
Sloan Digital Sky Survey (SDSS). We have shifted these by 
ASSFR = +0.6 to account for differences in the normalisa- 
tion of the SFR and in the dust correction. 



7 THE SPECIFIC STAR FORMATION RATE 
IN DIFFERENT ENVIRONMENTS 

To study the influence of the environment on the SSFRs of 
galaxies we investigate its behaviour for field galaxies and 
for galaxies in groups. Group membership in the if-selected 
MUNICS catalogue is assigned according to a modified ver- 
sion of the friends-of-friends algorithm, specifi cally designed 
to co pe with photometric redshift datasets l|Botzler et al.l 
2004). In brief, the algorithm works in redshift slices, and 
corresponding structures in adjacent slices are combined a 
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Figure 9. Distribution of rest-frame B — K colours as a function of selection band, luminosity, and redshift. Upper left-hand panel: 
All galaxies with ^ z 1 in the if-band selected (solid line), /-band selected (dot-dashed line), ij-band selected (dashed line) and 
B-band selected (dotted line) MUNICS samples, divided into three different _R"-band luminosity classes. All other panels: Sub-divided 
into galaxies with ^ z ^ 0.5 (dashed line) and 0.5 ^ z ^ 1.0 (dotted line) for MUNICS-B {upper right-hand panel), MUNICS-R {lower 
left-hand panel), and MUNICS_K {lower right-hand panel). In each plot, the distribution of Poisson errors is also shown below the object 
histograms. Furthermore, typical B — K error bars are also indicated in every plot. 



posteriori. Additionally to the traditional linking parame- 
ters in projected distance (Dl) and velocity space (Vl), a 
redshift pre-selection is performed to exclude objects with 
large photometric redshift errors. The linking criteria used 
for the MUNICS sample are D L = 0.125 h' 1 Mpc and 



Vl = 1000 km s 1 ; a minimum n umber of 3 objects is re- 
quired to form a group or cluster (|Botzler et al. 1 120071) . Ex- 
pected completeness fractions and contamination rates are 
roughly 90 per cent and 40 per cent, respectively. 

The resulting structure catalogue is presented in lBotzlerl 
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Figure 10. SSFR versus stellar mass for galaxies in MUNICS_K (symbols). Objects are coloured according to the age of the composite 
stellar-population synthesis model fit to their photometry, ranging from 0.01 Gyr (purple) to 8 Gyr (red). The solid and dashed black 
lines correspond to SFRs of 1 Mq yr — 1 and 5 Mq yr _1 , respectively. The dot-dashed line is the SSFR required to double a galaxy's mass 
between each redshift epoch and today, assuming constant SFR. The corresponding look-back time as well as an estimate of the typical 
error bars are indicated in each panel. Selection effects in this diagram are visualised for MUNICS_K (red line), MUNICS-I (orange line), 
and MUNICS.B (blue line). 



i|2004 ) and iBotzler et all (120071 1 and comprises 162 struc- 
tures (mostly groups) containing 890 galaxies in total. This 
group sample has already been used to study the integrated 
SSF R of groups in comp arison with field galaxies and clus- 
ters l|Feulner et al.ll2006h . 

A comparison of field galaxies and group members in 
the SSFR-stellar mass plane in two redshift intervals is pre- 
sented in Figure 1121 In contrast to the field population, 
group galaxies tend to populate the region of passive galaxies 
below the doubling line only. This, of course, is in agreement 
with previous findings that galaxies in higher density envi- 
ronments tend to have lower star-formation activity. Inter- 
estingly, however, passive galaxies are by no means limited 
to the group environment. This could be partially explained 
by group s election effects, but is in agreement with recent 
findings by iGerke et al.l (|2007l ) who identified a significant 
population of red galaxies in the field, with morphologies 
characteristic for passive galaxies. 

This behaviour is also obvious from Figure [13] where we 
show the redshift evolution of the average SSFR of galaxies 
in two different mass intervals, comparing field and group 
galaxies. At z ~ 1 there seems to be no marked difference 



between field galaxies and galaxies in groups, with a ten- 
dency of the difference increasing with decreasing redshift 
in the sense that group galaxies are rather more inactive 
than field galaxies. Note that this difference is larger for less 
massive galaxies than for the most massive objects. This is 
in principal agreement with earlier findings on the depen- 
dence of SSFR on dens ity for galaxies in the local universe 
l|Kauffmann et al.ll2004h . 

The fact that we do not see a more pronounced dif- 
ference between field and groups could be at least partially 
attributed to two effects. First, incompleteness of the group 
catalogue as well as contamination with false structures is 
likely to reduce the signal. Secondly, trends of the photomet- 
ric redshift errors with SFR and redshift might affect the 
analysis, especially in the higher redshift bins. Photometric 
redshift errors in the MUNICS catalogue increase with red- 
shift (from ~ 0.08 at z = 0.5 to ~ 0.12 at z = 1) and with 
SFR (they are typically 20 per cent larger for star-forming 
galaxies as compared to passive galaxies) . In reality, the dif- 
ference in SSFR between groups and the field may be more 
significant. 
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Figure 12. SSFR versus stellar mass for galaxies in MUNICS_K (symbols) for field galaxies (upper panels) and galaxies in MUNICS 
groups (lower panels). See caption of Figure flOl for details. 
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8 EVOLUTION OF THE LUMINOSITY 
FUNCTION 

8.1 Computing the Galaxy Luminosity Function 

As for t he near-infrared lum inos ity function (LF) pre- 
sented in iFeulner et~ail (|2003h and IDrorv etafl (|2003l ). the 
LF is compute d using the non-parametric Vmsx formalism 
|Schmidti ri968V In brief, the Vmax formalism accounts for 
the fact that some fainter galaxies are not visible in the 
whole survey volume. In a volume-limited redshift survey, 
each galaxy i in a given redshift bin [ziowor, z upP er] con- 
tributes to the number density an amount inversely pro- 
portional to the volume Vi of the survey: 

Zupper 

Vi = j (1) 

z lowoi 

where dV = dfl r 2 dr is the co- moving volume element 
and dfi is the solid angle covered by the survey. However, 
due to the fact that we have to deal with magnitude-limited 
surveys, a faint galaxy may not be visible in the whole sur- 
vey volume. Assuming that in a survey with given limiting 
magnitude galaxy i can be seen out to redshift z ma x, we have 
to correct the volume factor by V™"* /Vi, where 

mm(z uppcr ,z max ) 

vr-= J d ^d, ( 2) 

z lowcr 

Obviously, we have z max ^ z uppcr for a galaxy which 
is bright enough to be seen in the whole volume in in- 
vestigation, and the correction factor is one. Otherwise, 
Zmax < Supper, and the volume is smaller than the volume 
corresponding to the redshift range in which we compute the 
LF. We have made sure that the effect of the volume cor- 
rection is of importance only in the faintest bin in absolute 
magnitude, and that even in this case the correction is at 
most a factor of 5. 

Additionally, the contribution of each galaxy i is 
weighted by the inverse of the detection probability 
P(mdct,i) °f the parent catalogue. The LF $(M) is then 
computed according to the formula 

MM) dM = V i — dM, (3) 

V ; ^ Vi P(m doM ) ' V ; 

i 

where the sum runs over all objects i in the redshift 
range for which we want to calculate the LF. Naturally the 
volume terms can be simplified to l/V™* . 

We include two effects in the total error budget of the 
luminosity function. Firstly, the limited number of objects in 
each magnitude bin produces statistical uncertainties. Sec- 
ondly, the errors in the photometric redshift estimates have 
some influence on the luminosity function. 

The statistical errors are derived u sing Poi s sonia n 
statistic s, follow ing t he metho d s desc ribed in lGehrelsl l|l986h . 
lEbelind (|2003h . and lEbelind (|2004T ) for approximation of 
Poissonian errors for small numbers. 

To investigate the influence of photometric redshift er- 
rors on the luminosity function, we perform Monte-Carlo 
simulations. In these simulations, we take the original in- 
put catalogue used for deriving the luminosity function, 
and assign to each object a redshift within the redshift 



error distribution given by the photometric-redshift algo- 
rithm. In principle, there are two ways of doing this. Either 
one uses the redshift probability distribution of the best- 
fitting SED, or the sum of the probability distributions of 
all SEDs (the total probability distribution) . One might ex- 
pect the total probability distribution to have the advan- 
tage that it accounts better for systematic uncertainties in 
the SED fitting. However, we have carried out careful tests 
which show that the errors derived from the total distri- 
bution are not significantly different from the ones for the 
best distribution. This is due to the fact that the probabil- 
ity distribution for the best-fitting SED usually dominates 
the total distribution. Since the total distribution also takes 
more Monte- Carlo realisations to converge, we decided to 
use the probability distribution of the best-fitting SED for 
our simulations. The median of the redshift error assigned 
by the photometric redshift algorithm to galaxies in MU- 
NICS_R is ((j) mcd ian — 0.08, in good agreement with the 
measured deviation of photometric and spectroscopic red- 
shifts of Sz/(1 + z) ~ 0.057. Hence these error distributions 
are a good representation of the true error distribution. 

The errors are then computed as follows. The Poisson 
error and the standard deviation around the mean derived 
from the Monte-Carlo simulations are summed quadrati- 
cally. In addition, any difference between the measured value 
for the LF in a magnitude bin and the mean from the Monte- 
Carlo simulations is considered a measure of systematic er- 
rors. This is also added quadratically, but only in one direc- 
tion, i.e. if the Monte-Carlo mean is higher than the mea- 
sured value the upper error bar is enlarged, but not the lower 
one. All our \ 2 fitting routines, both for the Schechter pa- 
rameters and for the LF evolution, can handle asymmetric 
errors. 

8.2 i?-Band Selected Luminosity Functions 

In this Section we present luminosity functions from the R- 
selected MUNICS sample in the redshift intervals [0.4,0.6], 
[0.6,0.9], [0.9,1.5], and [1.5,2.5]. The luminosit y function 
resul ts together with a Schechter approximation l|Schechterl 
Il976h and a local comparison function are presented in Fig- 
ure 03] for the B band and in Figure [TS] for the R band, 
respectively. We also compare our r esults to the LF derived 
in the FDF in lGabasch et ail (|2004l2006l ). The agreement is 
very good in general, although we see a slight excess of very 
bright galaxies in MUNICS, resulting in somewhat larger 
characteristic luminosities and, because of the degeneracy 
of the Schechter paramters <&* and M* , lower characteristic 
densities. The excess of bright galaxies might be caused by 
three effects. First, quasars which cannot be easily discarded 
from our photometric catalo gue are expected to popula te the 
bright end of the LF (e.g. I Jahnke fc Wisotzkil [20031 ) . Sec- 
ondly, a slight excess of bright galaxies as compared to the 
Schechter function has bee n noted in spectroscopic surveys 
in th e local universe (e.g. iBlanton et ail 120031 ; IJones et ail 
2006). Finally, it may be at least partially attributed to 
photometric redsh ift errors s (see, e.g., the simulations in 
IDrorv et aill2003h . although in this case the error bars at 
the bright end of the LF seem to be smaller than the true 
errors which might point to an underestimation of the errors 
for individual galaxies, although they are in good agreement 
with the measured errors on average. In the highest redshift 



16 Georg Feulner et al. 



-1 
-2 
-3 
-4 
-5 
-6 



-1 
-2 
-3 
-4 
-5 
-6 



i ■ 1 ■ i ■ 1 ■ i ■ 1 

- MUNICSR 0.40<z<0.60 
MGC, Driver etal. (2005) 

- FDF, Gabasoh et al. (2004) 




* =0.0042, M =-21 .39, <x=-1 .25. 



-26 -24 -22 -20 -18 -16 -14 
Mb [mag] 



i ■ 1 ■ i ■ 1 ■ i ■ 1 

- MUNICS R 0.90 < z < 1.50 
MGC, Driver etal. (2005) 

- FDF, Gabasch et al. (2004) 




* =0.0023, M =-21 .74, <x=-1 .25. 



-26 -24 -22 -20 -18 -16 -14 
Mb [mag] 



o 

Cl 



o 



-1 
-2 
-3 
-4 
-5 
-6 



o 

Cl 



CO 




o 



-1 
-2 
-3 
-4 
-5 
-6 



i ■ 1 ■ i ■ 1 ■ i ■ 1 ■ i 

- MUNICS_R 0.60 < z < 0.90 
MGC, Driver etal. (2005) 

- FDF, Gabasch et al. (2004) 




* =0.0031 , M =-21 .47, oc=-1 .25. 



-26 -24 -22 -20 -18 -16 -14 
Mb [mag] 





I I I I I 
— MUNICS_R 1.50 <z< 2.50 






■■ MGC, Driver et al. (2005) 






- - FDF, Gabasch et al. (2004) 














/ / 
/ / 


t 






/ i 
I i 


<5*=0.0007, M*=-22.86, a= 


-1.25: 




. 1 / . 


I . , . I . , . I . , . I 





-26 -24 -22 -20 -18 -16 -14 
Mb [mag] 

0.50 (upper left-hand panel), z 



Figure 14. The B-band LF from the R-selected catalogue at redshifts z = 0.50 (upper left-hand panel), z = 0.75 (upper right-hand 
panel), z = 1.20 (lower left-hand panel) and z = 2.00 (lower right-hand pane l). One can clearly see the effect of brightening and 
decreasi ng number dens i ty wi th increasing redshift z. The local B-band LF from iDriver et ah (2005) is shown as dotted line, the FDF 
LF from lGabasch et al.l [120041) in the redshift bins [0.45,0.81], [0.45,0.81], [0.81,1.11] and [2.15,2.91] as dashed line. 



Table 3. Schechter parameters and M* and a for the LFs in the B and R band derived from MUNICS-R. 



Band 


<*> 




M* 


a 




min (x? cd ) 






[10" 3 Mpc" 3 ] 


[mag] 




[mag] 


B 


0.50 


4.20 ±0.60 


-21.39 ± 0.18 


-1.25 


-18.0 


0.37 




0.75 


3.12 ± 1.08 


-21.47 ±0.25 


-1.25 


-20.5 


1.11 




1.20 


2.32 ± 1.28 


-21.74 ±0.29 


-1.25 


-20.5 


0.28 




2.00 


0.68 ±0.22 


-22.86 ± 0.22 


-1.25 


-22.0 


0.41 


R 


0.50 


2.96 ±0.36 


-22.64 ±0.18 


-1.33 


-18.0 


0.37 




0.75 


3.12 ±0.62 


-22.54 ±0.18 


-1.33 


-20.0 


0.58 




1.20 


1.16 ±0.60 


-23.18 ± 0.30 


-1.33 


-21.5 


0.28 




2.00 


0.72 ±0.30 


-23.66 ± 0.26 


-1.33 


-23.0 


0.22 
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Figure 15. The i?-band LF from the i?-selected catalogue at redshifts z = 0.50 (upper left-hand panel), z = 0.75 (upper right-hand 
panel), z = 1.20 (lower left-hand panel) and z = 2.00 (lower right-hand panel). One can clea r ly see the effect of brightening and 
decreasing number density with increasing redshif t z. The local r'-band LF from iBlanton et alj l|2003l . transformed to the R band) 
is shown as dotted line, the FDF r-band LF from Gabasch ct al. (2006, transformed to the R band) in the redshift bins [0.45,0.85], 
[0.45,0.85], [0.85,1.31] and [1.91,2.61] as dashed line. 



bin, the MUNICS data do not sample the knee of the LF 
very well, resulting in rather low values for $*. 



8.3 Luminosity Function Evolution 

To estimate the rate of evolution of the Schechter parameters 
$* and M* with redshift, we define evolution parameters a 
and b as follows: 



We summarise the parameters M* and a of the 
Schechter fits in Table [3] and show the corresponding 
Schechter contours in Figure 1161 Note that we have kept 
the faint-end slope a fixed to the value —1.25 (B band) 
and —1.33 (R band) during the fitting process. Thes e are 
the values derived for the FDF in lGabasch et al.l (|2004h and 
iGabasch et ail (2006), respectively, which are in good agree- 
ment with the MUNICS LFs in the lowest redshift bin. At 
higher redshifts, the MUNICS data are not sufficiently deep 
to constrain the faint-end slope. 



M* (z) = M* (0) + oln(l + z) , 

$* (z) = $* (0) (1 + zf , and (4) 
a(z) — a (0) = a. 

Our parametrisation i s identical to the one chosen in 
IGabasch et "all (|2004 1 20061) and equivalent to the form $ oc 
(1 + z) p , L oc (1 + z)® sometimes found in the literature 
(especially in the context of radio-source evolution). The 
parameters translate as follows: 
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Figure 16. Schechter contours in the $*-M* plane for the ZJ-band LF (left-hand panel) and the ij-band LF (right-hand panel). 
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Figure 17. The redshift evolution of the luminosity function in MUNICS-R for the B-band LF (upper panels) and the ij-band LF (lower 
panels). In each case we show <&* versus z (left-hand panel) and M* versus z (right-hand panel) together with the evolutionary model 
fit described in the text (solid line). Note that this is not the fit to the Schechter parameters given in Table [3] and represented by the 
filled squares, but a simultaneous fit to t he LF values in all ma gnitude and re dshift bins (see text f or details). The dashed lines denote 
the evolution as derived from the FDF in lGabasch et aL | ||2004| , B band) and lGabasch et al.l ||2006| . r' band, Case 3). 
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and b — P. 
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Note th at this evolutionary m odel is different from th e 
one used in iFeulner et ail l|2003j ) and iDrorv et all ||200&D . 
where we used a parametrisation linear in redshift, $* (z) = 
<E>*(0) (1 + cz) and M*{z) = M*(0)+dz with the evolution- 
ary parameters c and d. In this case the parameters trans- 
late in the following way: a = dzj ln(l T z)~<i(z<l) and 
b — ln(l + cz)/ ln(l + z) ~ c (z <C 1). The advantage of the 
model used in this paper is that it is a good representation 
of the redshift evolution of $* and M* even at higher red- 



shift l|Gabasch et al.l 120041 ) . At low redshifts, however, the 
evolutionary models are equivalent. 

Since there is still some debate about the local galaxy 
LF (see e.g. the discussion in iDriver et alH2005l ), we decided 
to derive the evolutionary parameters from our MUNICS_R 
measurements alone. This means that we have to obtain the 
best-fitting values for a, b, $*(0), and M*(0) by minimising 
the four- dimensional x 2 distribution. 

The resulting values for the evolutionary parameters a, 
b, $*(0) and M*(0) can be found in Table Q] We show the 
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Figure 18. Estimates for the evolutionary parameters b and a 
(1(t and 2(r contours) for the B- and ij-band LFs computed on 
MUNICS-R. The solid lines show constant <E>* and M* , respec- 
tively, the dashed line indicates constant luminosity density, i.e. 
b = islga ~ o.921a. 




'pi = const. 



_i i i i i i i i i i_ 



M = const. 



_i i i i i i_ 



error contours of the evolutionary parameters a and b for the 
_R-selected MUNICS catalogue and the LF in the rest-frame 
B and R filter in Figure [TH] These contours were derived by 
projecting the four- dimensional \ 2 distribution to the a-b 
plane, i.e. for given a and b we used the values of M*(0) and 
$*(()), which minimise x 2 ( a i^)- 

We find 2a evidence for evolution in $* and M* for both 
filters. Furthermore, the evolution is stronger in B than in 
R, in agreement with the findings of iGabasch et all l|2004l . 
2006). In Figure [TBI we also show the line of constant lumi- 
nosity density pL- The evolution of the LF parameters seems 
to follow this relation rather closely, with la evidence for a 
slight decrease of the luminosity density with redshift.. 

The same evol utionary trend for the LF ca n be seen out 
to higher redshifts: [Gabasch et all (12004120061) use the deep 
/-band selected FORS Deep Field (FDF; lHeidt et al.ll2003l ) 
to trace the evolution of the LF from the ultraviolet to the 
near-infrared out to redshifts z ~ 5 finding similar results 
to the ones presented here. That dataset is complimentary 
to our catalogue, since the FDF lacks the large local volume 
of MUNICS which allows us to study the lower redshift part 
of the evolution with high statistical accuracy. 

In Figure [17] we show the resulting evolutionary model 
for $* and M* as a function of redshift, together with the 
Schechter parameters derived for the LF in the four redshift 
bins, local comparison values from the Millennium Galaxy 
Catalogue (|Driver et al.ll2005l ) for the B band and from the 
Sloan Digital Sky Survey (SDSS) presented in lBlanton et al.l 
ll2003|) for the R band, and the ev olution as derived from the 



9 SUMMARY AND CONCLUSIONS 

In this paper we presented results on the evolution of field 
galaxies drawn from galaxy catalogues selected in the B, R, 
I, and K filters in the context of the Munich Near-Infrared 
Cluster Survey (MUNICS). The first part of this paper is 
dedicated to a discussion of the construction and properties 
of the B, R, and /-selected MUNICS catalogues contain- 
ing ~ 9000, ~ 9000, and ~ 6000 galaxies, respectively. The 
catalogues reach 50% completeness limits for point sources 
of B ~ 24.5 mag, R ~ 23.5 mag, and / ~ 22.5 mag and 
cover an area of about 0.3 square degrees.. Photometric 
redshifts are derived for all galaxies with an accuracy of 
5z/(l + z) ~ 0.057, as demonstrated by comparing them 
with a sample of ~ 600 spectroscopic redshifts available for 
MUNICS galaxies. Star-galaxy classification is performed by 
comparing the \ 2 values of the spectral-energy-distribution 
fitting within the photometric redshift code. A comparison 
with spectroscopy clearly shows the reliability of this ap- 
proach. 

A study of rest- frame B — K colour distributions for the 
four different selection bands and different redshifts allows 
important conclusions about selection effects and the evo- 
lution of galaxy populations. First, we could show that the 
selection band influences the colour distributions only for 
objects with low A-band luminosities (comparatively low 
stellar masses) , while the distributions look remarkably sim- 
ilar for high-mass galaxies. Secondly, the colour distributions 
get wider and bluer with decreasing A luminosity, indicat- 
ing a higher average star-formation rate and a wider dis- 
tribution of star-formation rates for less massive galaxies. 
Thirdly, there is strong colour evolution with redshift for 
the redshift interval < z < 1 in the sense that with in- 
creasing redshift galaxies become bluer (i.e. have larger star- 
formation activity). This can hardly be seen for the most 
massive galaxies, but becomes more and more so for the 
lowest mass galaxies. Also, the trend becomes stronger go- 
ing from A-band selection to /3-band selection. This means 
that the increase in star formation rate from redshift zero 
to one is largely driven by lower mass galaxies, and that 
the most massive galaxies have assembled the bulk of their 
stellar mass before redshift unity, in agreement with our 
results on the specific star formation rat e (the star forma- 
tion rate per unit stel lar mass; see, e.g., Bauer et al. I l2005l ; 



FDF lGabasch et all (|2004l. 120061 ). 



iFeulner et"aTll2005l a.b ; lJuneau et al. 2005) 

We investigate the influence of selection band and en- 
vironment on the specific star formation rate (SSFR). We 
find that A-band selection indeed comes close to selection in 
stellar mass, while /3-band selection purely selected galaxies 
in star formation rate. As a second order effect, the depth of 
the other filters of the survey influence the selection bound- 
ary through the model fitting to the objects' photometry in- 
volved in deriving stellar masses. We use a galaxy group cat- 
alogue constructed on the A-band selected MUNICS sample 
to study possible differences of the SSFR between the field 
and the group environment, finding a marginally lower av- 
erage SSFR in groups as compared to the field, especially at 
lower redshifts. 

The field-galaxy luminosity function in the B and 
R band as derived from the /?-selected MUNICS cata- 
logue changes with increasing redshift in the sense that 
the characteristic luminosity increases but the number den- 
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Table 4. Parameters a, 6, M*(0), and 5>*(0) for the evolution of the LFs in the B and R band derived from MUNICS.R. 



Filter 


b 


a 


3>*(0) [10" 3 Mpc~ 3 ] 


Af*(0) [mag] 


ct 


min (x? cd ) 


B 


-2.47 ±0.58 


-2.07 ±0.38 


13.3 ±4.2 


-20.29 ±0.27 


— 1.25 


0.37 


R 


-1.84 ±0.61 


-1.43 ±0.43 


7.0 ±2.1 


-21.89 ±0.27 


— 1.33 


0.33 



sity decreases. This effect is smaller at longer rest-frame 
wavelengths and gets more pronounced at shorter wave- 
length s. This evolutionary trend continues to higher red- 
shifts (|Gabasch et al. 1 12004 l200rj ). Parametrising the red- 
shift evolution of the Schechter parameters as M*(z) = 
M*(0) + a ln(l + z) and $*(z) = $*(0)(1 + z) b we find 
evolutionary parameters a ~ —2.1 and b ~ —2.5 for the £? 
band, and a ~ — 1.4 and b ~ — 1.8 for the i? band. 
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APPENDIX A: PHOTOMETRIC REDSHIFTS 
IN DIFFERENT SAMPLES 

It is an important consistency check to compare the photo- 
metric redshifts of an object present in two different selec- 
tion catalogues. Note that the two redshifts do not have to 
be the same in the case of MUNICS, since small differences 
in object centring and differences in the object's shape in 
the different detection images lead to slightly different pho- 
tometry. Direct comparison of magnitudes in the same filter, 
but derived in the different selection catalogues sho w excel- 
lent agreement (comparison plots are published in iFeulnerl 
l2004h . 

Due to this small differences in the objects' photometry, 
the result for a comparison of the photometric redshifts in 
different catalogues is a scatter around the one-to-one rela- 
tion, as shown in Figure [1] However, the general agreement 
is extremely good, as can be seen from the histograms of 
redshift differences also shown in the Figure. The deviation 
from a zero redshift difference is smaller than 0.01 in z, and 
the a of the best-fitting Gaussian is smaller than 0.02, which 
is the redshift binning size of our photometric redshift algo- 
rithm. 
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Figure 1. Left-hand panels: Comparison of photometric redshifts from the /-selected MUNICS catalogue to those from the i^-selected 
MUNICS sample. Middle panels: The same for the R and the /-selected catalogue. Right-hand panels: The same for the B and the 
R-selected catalogue. The lower panels show histograms of the redshift differences together with a Gaussian fit showing the quality of 
the photometric redshifts. 



